Probing vortices in 4 He nanodroplets 
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We present static and dynamical properties of linear vortices in 4 He droplets obtained from Den- 
sity Functional calculations. By comparing the adsorption properties of different atomic impurities 
embedded in pure droplets and in droplets where a quantized vortex has been created, we sug- 
gest that Ca atoms should be the dopant of choice to detect vortices by means of spectroscopic 
experiments. 
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The unique environment realized in liquid He clus- 
ters has opened up in recent years new opportunities for 
atomic/molecular spectroscopy to probe superfluid phe- 
nomena on the atomic scale [J, |2j • Helium droplets rep- 
resent ideal nano-scale cryostats for a variety of funda- 
mental experiments on liquid 4 He, including the study of 
quantized vortices Q. Vortices, while energetically un- 
favorable Q, can potentially be stabilized by atomic or 
molecular impurities 

During the free jet expansion experiments described in 
Refs. |3j , it is plausible that quantized vortices may 
be created in some metastable state, long-lived enough 
to be detected. However, the question of whether 4 He 
droplets can sustain vortices is still not resolved, and all 
the high resolution spectra of embedded molecules can 
apparently be explained without invoking their presence. 
Yet, it is expected that in the near future they could 
be created by some extension of the present experimen- 
tal techniques. This calls for identifying signatures that 
might reveal vortical states in helium droplets. A pos- 
sible experiment to detect their presence has been de- 
scribed by Close et al. 0. They have suggested that 
alkali atoms, that normally reside in a 'dimple' on the 
surface of 4 He clusters 0, 0: ma Y be drawn, when 
a vortex is present, inside the cluster along the vortex 
core. Spectroscopic experiments on the dopant atoms 
could thus provide evidence of their existence, since the 
line broadenings and shifts would be different in the two 
cases. We show in the following that alkali atoms are 
actually not suited for such an experiment, but rather 
alkaline earth (Ca) atoms may serve as probes to detect 
vortices. 

Density Functional (DF) methods have become in- 
creasingly popular in recent years as a useful computa- 
tional tool to study the properties of classical and quan- 
tum inhomogeneous fluids, especially for large systems, 
for which they provide a good compromise between ac- 
curacy and computational cost. In particular, a quite 
accurate description of the properties of inhomogeneous 
liquid 4 He at zero temperature has been obtained within 



a DF ap pro ach by using the energy functional proposed 
in Ref. ^3] and later improved in Ref. [ijj. This later 
DF, which has been successfully used over recent years 
to study a variety of 4 He systems like clusters and films, 
is the one we use in the present work. 

The minimization of the energy functional with respect 
to density variations, subject to the constraint of a given 
number of 4 He atoms TV, leads to the equilibrium parti- 
cle density profile p(r), thus allowing to study the static 
properties of the 4 He system. When dynamical proper- 
ties are studied (as described in the following), we use 
the Time-Dependent DF (TDDF) method developed in 
Ref. 0], which allows to obtain both the 4 He particle 
density p(r, t) and the velocity field v(r, t). Briefly, in the 
static (dynamic) case, one has to solve a stationary (time- 
dependent) non-linear Schrodinger-like equation for an 
'order parameter' ^(r) (^(r, t)), where the Hamiltonian 
operator is given by H — -^V 2 + U[p,v]. The effec- 
tive potential U is defined as the variational derivative of 
the energv_ functional, and its explicit expression is given 
in Ref. [IJ. From the knowledge of \& = 0e* e one can 
get the density p(r, t) — <fi 2 and the fluid velocity field 
v(r, t) = — VO. To model the interaction of liquid 4 He 
with foreign impurities we use suitable He-impurity pair 
interaction potentials, which will be described later on. 

We work in 3D cartesian coordinates, and adopt the 
following procedure to generate a vortex in the cluster in 
the most unbiased way. We consider a cluster in a rotat- 
ing frame of reference with constant angular velocity lo z 
around the z-axis . The Hamiltonian density H then 
acquires an additional term —lo z L z , L z being the angular 
momentum component along the z-axis. We minimize $ 
for this constrained Hamiltonian, imposing "J to be or- 
thogonal, during the minimization, to Wo = v/ Peg ( r ) de- 
scribing the minimum energy state of the vortex-free clus- 
ter; we have applied the method to a TV = 300 droplet. 
To generate a vortex line, uj z must be larger than a crit- 
ical value -unknown in advance- Q c — AE/(NH) [T^| . 
where AE is the energy cost to create a vortex (which 
in the present case is about 70 K, see Table U and hence 
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51 c ~ 3 x 10 10 s^ 1 ), but not so large that one could 
generate a vortex array [T^- The particle density corre- 
sponding to this vortical configuration is shown in Fig. 
1(a). We have calculated the circulation of the velocity 
field along a path enclosing the vortex core, and have ex- 
actly found the value Nh/m appropriated for a quantized 
vortex with (L z ) = Nh. 

Note that since the vortex is quantized, the vortical 
state is an eigenstate of the angular momentum along the 
rotation axis, L z . This means that our density profile is 
the same as that one would obtain by using the Feynman- 
Onsager ansatz, i.e. by adding to the energy functional 
an extra centrifugal term associated with an order pa- 
rameter of the form ^/pe 1 ^ (<& being the azimuthal angle), 
and finding the density profile by solving an equation in 
the real quantity yj p(r). This is the procedure used in 
the DF calculations of Ref. [B| to generate quantized vor- 
tex structures in helium drops -and also in Bose-Einstein 
condensates of trapped gases ^ij-. Instead, we have not 
assumed a priori a quantized value for the total angu- 
lar momentum, but rather we generate a fully quantized 
vortex state starting from a pure cluster. 

To use atomic impurities as probes of the presence of 
vortices in 4 He drops, ideally one would like to have an 
atom that is barely stable on the surface of a pure drop, 
and becomes solvated in its interior in the presence of a 
vortex. The question of solvation vs. surface location for 
an impurity atom in liquid 4 He can be addressed in an 
approximate way within the model of Ref. |l5j | where, 
based on calculations of the energetics of impurities in- 
teracting with liquid 4 He, a simple criterion is proposed 
to decide whether surface or solvated states are favored. 
An adimensional parameter is defined in terms of the 
impurity-He potential well depth e and the minimum po- 
sition r m , A = p e r m /(2 1 ' 6 (r), where p and a are the 
bulk liquid density and surface tension of 4 He, respec- 
tively. The criterion for solvation reads A > 1.9 for the 
existence of solvated states [l5|- One thus needs an im- 
purity with A ~ 2, and such that its most stable state is 
on the drop surface. 

Alkali atoms are known to have their stable state on 
the surface of liquid 4 He [aH, and they lie in the low 
A regime (A ~ 0.6 — 0.9) [lq|. Accordingly, for alkalis 
a surface state should always be preferred, even in the 
presence of a vortex line. We have verified this point 
considering Na and Rb, as representative of light and 
heavy alkali, respectively. The alkali-Hc interaction is 
of the form proposed by Patil 0. We have compared 
the stable 'dimple' states of Na and Rb atoms on the 
surface of the N = 300 cluster hosting a vortex line, with 
those of the same impurity trapped in the vortex core, 
exactly at the cluster center. We have found that the 
latter are energetically unfavored with respect to surface 
states, see Table [I] It is worth to note that in the case 
of Na, our results compare well with the Path Integral 
Monte Carlo calculations of Ref. ^tJi where a binding 
energy of about ~ 7 K is found for this cluster. We also 
note that, unlike the case of strongly bound impurities 



to 4 He clusters which have their stable state inside 
the cluster, and for which there exists a critical cluster 
size below which the droplet+dopant+vortex complex is 
stable, the alkalis cannot stabilize the vortex, whatever 
the droplet size is 0- 

There are other dopants, however, for which there is 
clear evidence of a surface state on liquid 4 He, i.e. al- 
kaline earth atoms. Absorption spectra of alkaline earth 
atoms (Ca, Ba and Sr) attached to 4 He clusters clearly 
support an outside location of Ca and Sr [Tsj , and proba- 
bly also of Ba 0] . To describe the He- impurity interac- 
tion we employ an accurate ab-initio He-Ca pair potential 
[20| used to study 4 HeAr+Ca droplets up to TV — 75 by 
Diffusion Monte Carlo techniques |2l|. For such a po- 
tential, A ~ 2.2, which apparently indicates a solvated 
stable state. However, for those cases where A is close to 
the solvation threshold, consideration of the shape of the 
potential energy surface, as well as the well depth and 
equilibrium internuclear distance, seems warranted |22j. 

The stable state of a Ca atom in a 4 He3oo vortex-free 
cluster is shown in Fig. 1(b). Note that, in qualitative 
agreement with the experimental evidence, the 'dimple' 
appears to be much more pronounced than in the case of 
alkalis reflecting the stronger He-atom interaction. 
However, in the presence of a vortex, the stable state 
is in the center of the cluster, as depicted in Fig. 1(c). 
The surface state in this case is unstable and, as the 
minimization proceeds, a Ca atom initially placed on the 
surface near the vortex core, is gradually drawn towards 
it and then sucked inside, eventually reaching the stable 
state in the center of the drop. The value of the angular 
momentum for the converged 4 He configuration is again 
(L z ) =NH. 

The response of the impurity atom to the different 4 He 
environments shown in Fig. 1 might be determined with 
spectroscopic measurements, allowing to detect the pres- 
ence of vortices: the observed linewidths and shifts of 
the excitation/emission spectra in the two cases shown 
in Fig. 1 should be very different, reflecting the 'bubble' 
environment in one case [Fig. 1(c)], and a more open 
environment in the other case [Fig. 1(b)]. Moreover, in 
the case of Fig. 1(b), bound- unbound transitions should 
be observed with a significant probability, thus implying 
a strong asymmetry in the observed spectra. 

In the above picture one is assuming that the vortex 
is long-lived enough to allow a Ca atom, picked up ran- 
domly by the cluster, to diffuse close to the top of the 
vortex core and then to be drawn inside. We have no 
direct proof of the stability of the cluster+vortex com- 
plex on experimental time scales. However, we have indi- 
cations that the cluster+vortex+dopant complex should 
be stable at least on the nanosecond time-scale. This 
conclusion comes from very long computing time simu- 
lations, using the TDDF method 0, to study the dy- 
namics of the clustcr+vortex+impurity complex. Dur- 
ing these simulations, the impurity atoms were allowed 
to oscillate inside the 'bubble' in the cluster center [see 
Fig. 1 (c)] , and the vortex line was always found to be sta- 
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ble, without showing any tendency to shrink, bend or 
migrate towards the surface of the cluster. 

The solvated state for Ca in a vortex-free cluster is 
a stationary but unstable configuration against any dis- 
placement of the atom off the cluster center, only a few 
K in energy above the stable, 'dimple' state (see Table 
HJ . This is a consequence of the borderline value of A for 
this impurity, and implies that in a real experiment, a 
fraction of Ca atoms might be trapped inside the clus- 
ters for fairly long times, even in the absence of vortices. 
For these atoms, the spectroscopic signals would be simi- 
lar to those coming from Ca atoms trapped in the vortex 
core, making it difficult to discriminate between the two 
cases. A line shape calculation using as an input the 
4 He density profiles around the impurity might help to 
distinguish between solvated states of Ca with and with- 
out vortex. Since the extension of the method of Ref. |(j 
-which is applied there to the simpler case of the mono- 
electronic alkali atoms- is rather involved for two-electron 
systems, we have not carried out such a calculation. We 
instead suggest additional measurements which may help 
to discriminate between the states shown in Fig. 1. 

It appears from our calculations that the energy of 
the impurity-cluster system is rather insensitive to the 
location of Ca along the vortex core, once the atom is 
embedded in it. Consequently, vibrational modes of the 
impurity along the vortex line are expected to be soft. 
We have confirmed this by TDDF calculations, applying 
to the Ca atom a small initial momentum in a given di- 
rection -radially, towards the surface of the cluster for 
the 'dimple' state of Fig. 1(b), and along or perpendic- 
ular to the vortex line in the case shown in Fig. 1(c)-. 
We then let the impurity evolve in time, allowing for the 
4 He environment to dynamically follow the atom motion 
while the impurity oscillates around its equilibrium po- 
sition. This is done in practice by numerically solving 
by means of a discrete Verlet algorithm, as usually done 
in Molecular Dynamics calculations, Newton's equation 
of motion for the Ca atom under the force due to the 
surrounding 4 He liquid 

= -V R p(r,t)V He -ca(R(t) - r)drj 

In this expression Vne-Ca is the pair potential describ- 
ing the He-impurity interaction, and the density p{v, t) 
is updated at each time step according to the TDDF 
scheme for 4 He ^3 ■ From the positions of the Ca atom 
as a function of time, relative to the center-of-mass of the 
Ca-He droplet system, different frequencies characteriz- 
ing the impurity dynamics can be found from a Fourier 



analysis of the calculated time series. 

We report in Fig. 2 the calculated vibrational spectra. 
The intensities are in arbitrary units, and normalized so 
that the higher peak in each spectrum has unit height 
[23| . It appears that the oscillation of Ca along the vor- 
tex core is indeed characterized by a single low frequency 
mode, as compared with the more fragmented spectrum 
for vibrations perpendicular to the vortex core. The pres- 
ence of the soft mode is a signature of solvation of a Ca 
atom inside the vortex. Indeed, such a mode should be 
severely damped, or even absent, for a solvated Ca atom 
in a vortex-free cluster, since this configuration is unsta- 
ble. We also show for comparison the vibration spectrum 
of a Ca atom in the 'dimple' state on the surface of a clus- 
ter without vortex. The peak just below 1 K is due to 
the 'dipolar' vibration of the impurity inside the semi- 
spherical(spherical) cavity in which it is trapped in the 
'dimple' ('bubble') state. Additional peaks appear in the 
'dimple' spectra because of the coupling of the Ca motion 
with the surface modes of the 4 He nanodroplet, which 
have similar frequencies (for instance, the lowest energy, 
I = 2 quadrupolar mode of a pure 4 He3oo droplet occurs 
at ~ 0.6 K 12]). All these modes lie in the microwave fre- 
quency regime and there are experimental ideas to mea- 
sure the corresponding vibrational frequencies |24j . 

It is worth to see how these modes are coupled, which 
is particularly apparent when we displace the impurity 
perpendicular to the vortex core (dashed line in Fig. 2). 
To trigger this oscillation, we have given to the the impu- 
rity a kinetic energy of about 2 K, three times as much as 
in the other two cases. One may see that the spectrum 
displays one peak corresponding to the 'dipolar' mode 
discussed previously, and also softer modes of character- 
istics similar to those found in the other two cases. This 
coupling is possible because the time evolution is adia- 
batic. 

Finally, we would like to emphasize that the dynamical 
behavior of Hc-impurity systems depends in a sensitive 
way on the details of the He-atom pair interaction. This 
calls for improving the available interaction potentials to 
strengthen the scenario described here, or to help in find- 
ing other atomic/molecular impurities which may serve 
as probes of the presence of vortices in 4 He droplets. 
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Figure Captions 

Figure 1: From left to right, for a 4 Hc3oo cluster, 
atomic equidensity lines corresponding to: (a) the vor- 
tical state; (b) the deep 'dimple' stable state of Ca, and 
(c) the solvated stable state of Ca in the presence of a 
quantized vortex. Each box side has 45 A lenght. 

Figure 2: Calculated frequency spectra for the oscilla- 
tions of a Ca atom in a 4 He3oo cluster. Solid line: sol- 
vated state [see Fig. 1(c)], Ca moving along the vortex 
core; dashed line: solvated state, Ca moving perpendic- 
ular to the vortex core; dotted line: surface state [see 
Fig. 1(b)], Ca moving in the radial direction. 

TABLE I: Energies of several 4 He3oo configurations. In the 
left part of the Table, they are referred to the total energy of 
the pure, vortex-free cluster (-1384.5 K), whereas in the right 
part they arc referred to that of the clustcr+vortcx configura- 
tion (-1313.4 K). The configurations marked with an asterisk 
are unstable stationary configurations. 
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E (K) 


Rb (center)* 


+94.7 


Rb/vortex (center)* 


+71.7 


Rb (top) 


-9.1 


Rb/vortex (top) 


-12.5 


Na (center)* 


+64.3 


Na/vortex (center)* 


+45.1 


Na (top) 


-8.1 


Na/vortex (top) 


-11.3 


Ca (center)* 


-33.7 


Ca/vortex (center) 


-49.9 


Ca (top) 


-37.4 


Ca/vortex (top) 


unstable 


Vortex 


+71.1 







duration of our simulations, t max <■ 
[24] W. Ernst, private communication. 
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